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ABSTRACT 

Light curves from the Kepler Mission contain valuable information on the nature of the phenomena 
producing the transit-like signals. To assist in exploring the possibility that they are due to an 
astrophysical false positive, we describe a procedure (BLENDER) to model the photometry in terms of 
a "blend" rather than a planet orbiting a star. A blend may consist of a background or foreground 
eclipsing binary (or star-planet pair) whose eclipses are attenuated by the light of the candidate and 
possibly other stars within the photometric aperture. We apply BLENDER to the case of Kepler-9, a 
target harboring two previously confirmed Saturn-size planets (Kepler-9 b and Kepler-9 c) showing 
transit timing variations, and an additional shallower signal with a 1.59-day period suggesting the 
presence of a super-Earth-size planet. Using BLENDER together with constraints from other follow-up 
observations we are able to rule out all blends for the two deeper signals, and provide independent 
validation of their planetary nature. For the shallower signal we rule out a large fraction of the false 
positives that might mimic the transits. The false alarm rate for remaining blends depends in part 
(and inversely) on the unknown frequency of small-size planets. Based on several realistic estimates of 
this frequency we conclude with very high confidence that this small signal is due to a super-Earth-size 
planet (Kepler-9 d) in a multiple system, rather than a false positive. The radius is determined to be 
1.64+0 Rq, and current spectroscopic observations are as yet insufficient to establish its mass. 

Subject headings: binaries: eclipsing — planetary systems — stars: individual (Kepler-9, KOI-377) 
— stars: statistics 



1. INTRODUCTION 

The Kepler Mission, launched in March of 2009, was 
designed to address the important question of the fre- 
quency of Earth-size planets around Sun-like stars, and 
to characterize extrasolar transiting planets through a 
3.5-year program of ve ry precise photom etric monitor- 
ing of -156 000 stars (jKoch et al.ll2010D . Science re- 
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suits from the missio n have already begu n to appear 
(|Borucki et all I2010al |b1: ISteffen et all 1201 01 ) . As shown 
already by ground-based surveys for transiting planets, 
considerable effort is required to validate candidates de- 
tected photometrically. This is because false positives 
usually outnumber true planetary systems by a large fac- 
tor, which is about 10:1 for the most successful surveys 
from the ground, but is not yet well characterized for 
Kepler. The follow- up efforts by the Kep ler team have 
been summarized by iBatalha et al.1 ()2010f ) . 

Spectroscopy is often a crucial step in the vetting pro- 
cess, as it allows not only to measure the mass of a 
planet but also to examine any changes in the line pro- 
files (bis ector spans) that might indicate a fa lse posi- 
tive (see lOueloz et all 120011: iTorres et alll2005[ ). Some 
of the most challenging false positives to rule out in- 
clude chance alignments with a background eclipsing bi- 
nary ("blends"). However, for faint candidates (V > 14) 
high-resolution, high signal-to-noise ratio spectroscopy 
becomes prohibitively expensive in terms of telescope 
time. Even for brighter candidates, the reflex motion 
of the star due to an Earth-mass planet can sometimes 
be below the radial-velocity detection limit, making spec- 
troscopic confirmation very difficult or impossible. The 
question then becomes how to validate these candidates, 
particularly the ones with small planets that are precisely 
among the most interesting. 

A number of other tests have been developed that can 
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aid in understanding the nature of the candidate, and 
that rely on the long-term and nearly continuous pho- 
tometric monitoring of Kepler, as well as the very high 
astrometric precision achiev ed in determining the cen- 
troids of the stars (see also ISteffen et al1l2010l ). These 
tests include: i) verifying that alternating events have 
the same depth, which they may not if the signal is due 
to a background eclipsing binary; ii) checking for the 
presence of shallow secondary eclipses, which are com- 
mon in eclipsing binaries but are not expected for the 
smallest planets; Hi) checking for ellipsoidal variations, 
which could be another sign of a blend, iv) checking 
for changes in the centroid positions correlated with the 
brightness changes, which, if detected, might indicate a 
blend, or at the very least, a crowded aperture. This 
is a powerful diagnostic that is able to disprove many 
background blends. 

In addition to these tests, high-resolution imaging is an 
important tool to identify neighboring stars that might 
be eclipsing binaries with the potential to cause the 
transit-like signals. The photometric aperture of Kepler 
is large enough (typically many arc seconds across) that 
it usually includes other stars in addition to the candi- 
date, which increases the risk of such blends. In some 
cases, near-infrared observations with Warm Spitzer can 
allow one to reject the planet hypothesis if the transit 
depth at 3.6 /im or 4.5 /im is significantly different from 
that in the Kepler band. Such a signature might result 
from a blend with an eclipsing binary of a different spec- 
tral type than the candidate. 20 

Even with this extensive battery of tests it may still be 
difficult or impossible to provide validation for many of 
the most interesting planet candidates discovered by Ke- 
pler. For example, blend scenarios involving an eclipsing 
binary or an eclipsing star-planet pair physically asso- 
ciated with the candidate (hierarchical triple systems) 
and in a long-period orbit around their common cen- 
ter of mass would often be spatially unresolved from the 
ground. These configurations may also not be detectable 
spectroscopically, and would likewise not produce any 
measurable centroid motion. Therefore, it is imperative 
to take advantage of all the information available in vet- 
ting candidates. 

With this as our motivation, we describe here the use 
of the Kepler light curves themselves in a different way to 
help discriminate between true planetary transits and a 
large variety of possible blend scenarios, on a much more 
quantitative basis than simple back-of-the-envelope cal- 
culations could provide. The technique tests these sce- 
narios by directly modeling the light curves as blends, 
and has considerable predictive power that allows the 
expected properties of the various configurations to be 
tested against other information that may be available. 
Both hierarchical triples and background blends can be 
explored. A restricted application of this type of model- 
ing to Kepler has already been made for t he five multi- 
planet candidates announced recently by ISteffen et al.1 
(2010). For the present paper we have chosen to illus- 
trate the full potential of the method, which we refer to 

20 For Earth-size planets, the amplitude of the signal in the Ke- 
pler band is very small and possibly below the detection threshold 
for Spitzer. However, a blend with a late- type binary could pro- 
duce a much deeper eclipse at longer wavelengths that may be 
detectable in the near-infrared Spitzer bands. 



as BLENDER, by applying it to the unique case of Kepler 
Object of Interest 377 (KOI-377, henceforth Kepler-9). 
This is a multi-planet system reported and described in 
detail bv lHolman et al.l (|2010f ) . with three low-amplitude 
periodic signals in its light curve. We have selected this 
system for two main reasons. On the one hand, it rep- 
resents the first unambiguous detection of transit timing 
variations (TTVs) in an extrasolar planet, with a pat- 
tern of variation seen in two of its signals (Kepler-9 b and 
Kepler-9 c) that constitutes irrefutable evidence that the 
objects producing them are bona-fide planets. This offers 
an ideal opportunity to test BLENDER because their true 
nature is already known. On the other hand, the third 
signal (KOI-377. 03) 21 is small enough that it would cor- 
respond to a super-Earth, but validation of its planetary 
origin is not yet in hand. Should it be validated, Kepler- 
9 would become an even more remarkable laboratory for 
the study of the architecture of planetary systems in- 
volving small planets. Thus, exploring the wide range 
of possible blend configurations that might mimic this 
shallow signal is of the greatest interest for determining 
its true nature. 

Kepler is likely to find many other candidate transiting 
planets similar to KOI-377. 03, for which final validation 
by other means is not currently feasible, either because 
the expected radial-velocity signal is too small, or be- 
cause Doppler measurements are otherwise complicated 
due to the star being chromospherically active, rapidly 
rotating, or too faint. With the application to Kepler- 
9 we show that our light-curve modeling technique is a 
powerful tool for exploring astrophysical false positive 
scenarios that is complementary to other diagnostics, and 
should play an important role in the discovery of Earth- 
size planets around other Kepler targets. 

2. SIMULATING FALSE POSITIVES WITH BLENDER 

In general the detailed shape of a light curve display- 
ing transit-like events can be expected to contain use- 
ful constraints on possible blend scenarios that might be 
responsible for those signals. With photometry of the 
quality of that provided by Kepler, those constraints can 
be quite strong, and may be used to exclude many blend 
configurations and provide support for the planetary hy- 
pothesis. It is thus highly desirable to take advantage 
of this information, particularly since it relies only on 
observations already in hand. 

The idea behind BLENDER is to compare the transit 
photometry of a candidate against synthetic light curves 
produced by an eclipsing binary that is included within 
the photometric aperture of Kepler, and is contaminating 
the light of the candidate. The usually deep eclipses of 
the binary are attenuated by the light of the candidate, 
and reduced in depth so that they appear transit-like. 
In principle there is an enormous range of possible bi- 
nary configurations that could mimic all of the features 
of true planetary transits, including not only their depth, 
but also the total duration and the length of the ingress 
and egress phases. Generally it is only with detailed 
modeling that these can be ruled out. Possible scenarios 
include not only background eclipsing binaries, but also 

21 The name of this candidate follows the convention of the Ke- 
pler Mission in which individual transiting planet candidates are 
designated with a numerical tag, and validated planets are given a 
Kepler number and letter designation as in Kepler-9 b. 
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hierarchical triples, i.e., an eclipsing binary physically as- 
sociated with the candidate in a wide orbit around their 
common center of mass. 

The basic procedure for simulatin g light curves wit h 
BLENDER was described in detail bv ITorres et al.1 (|2004f ). 
and further changes and enhancements are discussed be- 
low. Briefly, the brightness variations of an eclipsing 
binary are generated with the binary light-curve code 
EBOP ([Popper fc Etzell Il981h. bas ed on the Nelson- 
Davis-Etzel model ( [Nelson & Davisl Il972t lEtzell I1981D . 
and then diluted by the light of the candidate for com- 
parison with the Kepler observations. Effects such as 
limb darkening, gravity brightening, mutual reflection, 
and oblateness of the binary components are included. 
The objects composing the binary are referred to as the 
'secondary' and 'tertiary', and the candidate is the 'pri- 
mary'. The properties of each object needed to gener- 
ate the light curves ( brightness and size) are taken from 
model isochrones by iMarigo et al.l ([2008!) , parametrized 
in terms of their stellar mass. 22 For the primary star the 
appropriate isochrone is selected by using as constraints 
the effective temperature, surface gravity, and metallic- 
ity determined spectroscopically. We assign also a mass 
and a radius from this isochrone, although these char- 
acteristics are irrelevant for generating the model light 
curves. We then read off the intrinsic brightness of the 
star (absolute magnitude) in the Kepler passband, which 
is the only property needed by BLENDER. The brightness 
of the primary is held fixed throughout all simulations. 
The parameters of the binary components are allowed 
to vary freely over wide ranges in order to provide the 
best match to the Kepler photometry in a chi-square 
sense, subject only to the condition that the two stars lie 
on the same isochrone, as expected from coeval forma- 
tion. To read off their properties (absolute magnitude 
and size) we use the mass as an intermediate variable. 
The specific isochrone adopted for the binary depends 
on the configuration: for hierarchical triple scenarios we 
adopt the same age and chemical composition as the pri- 
mary, whereas for background binaries the isochrone can 
be different. The Kepler light curve itself does not pro- 
vide a useful constraint on the age or metallicity of the 
binary in the background case, so a typical choice is a 
model for solar metallicity and a representative age for 
the field such as 3 Gyr. For background binary scenar- 
ios the distance between the binary and the main star is 
parametrized for convenience in terms of the difference 
in distance modulus, A5. The inferred distance between 
the primary star and the observer will vary from blend to 
blend because we constrain the combined brightness of 
all components of the blend to match the measured ap- 
parent brightness of the target. BLENDER is also able to 
account for differential extinction between the primary 
and the binary, which can have a non-negligible effect in 
some cases given the relatively low Galactic latitude of 



the Kepler field. 

Early versions of BLENDER have been used occasion- 
ally in recent years to examine transiting planet can- 
didates from ground-based surve ys such as O GLE , 
TrES , and HATNet ( see, e .g.. ITorres et all 12004 
20051 iMandushev et al.l 120051 : 10 'Donovan et al.l 120061: 



Bakos et al.ll2007D . as well as from CoRoT ( Fressin et al.l 
20101) . These studies have exploited the predictive power 
of BLENDER to estimate further properties of the blend 
scenarios, by testing them against complementary infor- 
mation such as color indices, optical/near- infrared spec- 
troscopy, or near-infrared photometry from Spitzer. For 
the application to Kepler, several important modifica- 
tions have been made to BLENDER, including the follow- 
ing: i) the ability to generate light curves integrated 
over the 29.4-minutc effective duration of an observation 
when using long-cadence data. This changes the shape of 
the transits significantly, given the high precision of the 
Kepler photometr y and the relatively short timescale s 
of the events (see iGilliland et afll2lnol : iKippind [20101) : 
ii) de-trending of the original Kepler light curves with 
a 1-day running median to remove instrumental effects, 
and rejection of outliers; Hi) the use of model isochrones 
specific to the Kepler passband, kindly computed for us 
by L. Girardi. BLENDER can now also use proper limb- 
darkening coefficients for the same band, as opposed to 
an approximation to the Kepler passband such as the 
Johnson R filter, which is considerably narrower. Ke- 
pler l imb-d arkening coefficients have been computed by 
ISingl (|2010f ) and also A. Prsa; 23 iv) extension to any 
optical or near-infrared passband. In particular, for any 
scenario explored with BLENDER, light curves can be com- 
puted at other wavelengths such as the 3.6 /im and 4.5 /mi 
passbands of the IRAC instrument on Warm Spitzer, in 
order to further test the blend hypothesis. Addition- 
ally, BLENDER can predict the overall color of a blend in 
any pair of passbands, including the effects of differential 
reddening for background or foreground scenarios. These 
colors may then be compared with the measured colors of 
a target. Extinction at differ ent wavelengths is co mputed 
following the prescription bv lCardelli et al.l (|1989h : v) the 
ability to have the tertiary be a (dark) planet instead of 
a star, in which case the corresponding free parameter 
becomes the radius of the planet rather than the tertiary 
mass. The mass of the planet has little effect on the light 
curves in most cases, but can nevertheless be set to any 
value; 24 vi) the ability to include extra light from other 
stars that may be present in the Kepler aperture, which 
further dilutes the intrinsic signatures from the eclipsing 
binary; vii) the ability to model systems with eccentric 
orbits. Eccentricity changes the orbital velocities during 
transit, and can therefore affect the size (mass) of stars 
that allow satisfactory fits to the light curve. 

When exploring blend scenarios involving hierarchical 
triple systems, the free parameters of the problem are the 



22 This particular set of isochrones was chosen because it reaches 
lower masses than other models (nominally 0.15 Mq, which we 
have extrapolated slightly for this application to 0.10 Mq, near 
the brown dwarf limit), and because a convenient web tool pro- 
vided by the authors allows easy interpolation in both age and 
metallicity (http://stev.oapd.inaf.it/cgi-bln/cmd). Addition- 
ally, isochrone magnitudes are available in a variety of passbands 
including the Kepler and Spitzer passbands, as well as Sloan and 
2MASS. 



23 http : //astro4 . ast . villanova. edu/aprsa/?q=node/8 

24 We note that this option of BLENDER implicitly allows to con- 
sider white dwarfs as tertiaries, as they are also Earth-size and 
contribute little light. Their mass is significantly larger than a 
planet's mass, however, which in close orbits can lead to distortions 
in the primary star causing ellipsoidal variation. Gravitational mi- 
crolensing may also occur in systems involving white dwarfs with 
long enough periods, and may well be detectable in the Kepler 
photometry. 
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mass of the secondary, the mass of the tertiary (or its ra- 
dius, if a planet), and the inclination angle. A fourth 
variable, the difference in distance modulus, is added 
for background blends. These quantities are typically 
stepped over wide ranges in a grid pattern to fully map 
the x 2 surface. For the application to Kepler-9 below, 
stellar masses are allowed to vary along the isochrones be- 
tween O.IMq and 1.4 M©, although at the larger values 
the observed duration of the transits is already difficult 
to match unless the events are highly grazing, in which 
case the shape would be very different. For planetary 
tertiaries the radii are allowed to be as large as 1.8 i?j up ; 
values higher than this have not been observed. 

3. APPLICATION TO KEPLER-9 

Kepler-9 (KIC 3323887, 2MASS 19021775+3824032) is 
a relatively faint star compared to typical ground-based 
transit hosts (Kepler magnitude Kp = 13.8), which was 
observed by the mission beginning in the first quarter 
of operations, and presents three distinct periodic sig- 
nals in its light curve. The two with the largest ampli- 
tudes have periods of 19.24 days (Kepler-9 b) and 38.91 
days (Kepler-9 c), and brightness decrements of 6.5 and 
6.0 mmag, respectively. The third signal (KOI-377.03) 
is much shallower (0.2 mmag), and repeats every 1.59 
days. The two longer periods are within 2.5% of being 
in a 2:1 ratio, and both objects display obvious TTVs 
that are anti-correlated, clearly indicating they are inter- 
acting gravitationally and therefor e orbit the same star , 
and are planetary in nature (see iHolman et ahl I2010D . 
The estimated radii are quite similar to that of Saturn, 
and the masses are somewhat smaller than Saturn, based 
on available radial- velocity measurements constrained by 
transit times and durations. The short-period signal has 
one of the smallest amplitudes detected by Kepler, and 
may well correspond to a third, super-Earth-size planet 
in the system, with a n estimated radius of only ~1.5 R® 
(jHolman et al.l 120101 ). However, because it shows no 
TTVs related to the other two planets (nor is expected 
to, on dynamical grounds), and is predicted to induce 
only a very small reflex velocity on the parent star that 
may be below detection for such a faint object, the true 
origin of this signal has not yet been established. 

In the absence of the crucial evidence of TTVs, each 
of the two largest signals — and indeed the third sig- 
nal as well — could in principle be due to a different 
blend. 25 Therefore, as an illustration of the application 
of BLENDER, we model the light curve of Kepler-9 at each 
period separately, as we would any candidate with a sin- 
gle period, and we account for possible blends at the 
other periods by incorporating extra dilution consistent 
with those other scenarios. The goal for the two largest 
signals is to demonstrate, as a sanity check, that BLENDER 
would be able to rule out blends in similar cases where 

25 Unlikely as it may seem to have three different blends oper- 
ating in the same system, the large photometric aperture, nearly 
uninterrupted monitoring, very high photometric precision, and 
long-term coverage of Kepler coupled with the large number of 
targets observed makes it more sensitive to picking up odd cases 
such as this, so they should not be completely ruled out. An exam- 
ple already exists among the fiv e multi-planet candidates recently 
reported bv lSteffen et al.l (2010[). in which one of the systems (KOI- 
191) presents three transit-like signals, and one of those signals (0.4 
mmag depth) has been shown to be due to a background eclipsing 
binary 2.6 mag fainter than the target, located 1.5 arcsec away. 



confirmation is lacking, which Kepler is expected to find 
in significant numbers. For the third signal of unknown 
nature, the application of BLENDER should provide valu- 
able evidence one way or the other. 

3.1. Stellar properties and photometry 

Kepler-9 is a solar type star. The spect roscopic prop- 
erties of the primary are adopted from IHolman et al.l 
(|2T)10t) : T cB = 5777±61K, log<? = 4.49±0.09, and [Fe/H] 
= +0.12 ± 0.04. With these para meters, a comparison 
with the stellar evolution models of iMarigo et al.l (|2008l ) 
yields a stellar mass of = 1.07+ 0.05 M , a radius of 
i?+ = 1.02 ± 0.05 Rq, and an age of about 1 Gyr, along 
with the absolute magnitude in the Kepler band. Only 
the latter is used by BLENDER, and is held fixed in our 
modeling. The distance to the star estimated from the 
same models is about 650 pc, ignoring extinction. Uncer- 
tainties in the brightness of the primary stemming from 
errors in T e g, logg, and [Fe/H] are small. For example, 
the error in log 5, which has the most direct influence on 
the intrinsic brightness, translates to an uncertainty of 
little more than 0.1 mag in the absolute magnitude. This 
has an insignificant impact on our results. 

The photometry used here consists of the long- 
cadence measurements gathered for Kepler-9 during Ke- 
pler quarters 1, 2, and 3, spanning 218 days, and was 
treated slightly differently than indicated earlier for a 
generic Kepler candidate because of the complications 
stemming from the TTVs. Using the binary FITS 
tables from MAST (Multimission Ar chive at STScI, 
http://archive.stsci.edu/kepler/), the "raw" aper- 
ture photometry for each quarter was first de-trended 
using a moving cubic polynomial fit robustly to out- 
of-transit data, with a sliding window of 999 minutes 
before and after each individual datapoint. This tech- 
nique removes long-term trends due to stellar activity or 
instrumental errors, but retains the properties of each 
transit light curve. Statistically significant outliers were 
removed. 

For the two long-period signals, simple folding will 
not create an accurate light curve because of the strong 
TTVs. Instead, we used a "shift-and-stack" technique, in 
which each transit event is displaced so that it is centered 
at "time" zer o usin g the measured transit times from 
IHolman et al.l (|2010| ). Along with the measurements in 
transit, nearly a full cycle of out-of-transit data were also 
shifted. Specifically, we shifted nearly 25% of an orbital 
cycle before the transit, and nearly 75% after the tran- 
sit. This preserves any curvature outside of eclipse, and 
in principle would also reproduce any secondary eclipses, 
both of which can provide useful constraints when model- 
ing the light curve with BLENDER. We note, however, that 
the strong TTVs of the transits would be accompanied 
by shifted secondary eclipses in a way that can only be 
predicted by full numerical integration. This shift-and- 
stack technique would not align secondary eclipses cor- 
rectly and thus their depth would need to be significant 
in each individual event to be noticed. There is no sign 
of secondary eclipses at these periods in the data at the 
10~ 4 level, as expected from the planetary nature of the 
objects, and thus the failure of the shift-and-stack tech- 
nique to correctly add up the secondary eclipses does not 
affect our results. After shifting, all the transit and out- 
of-transit data were "stacked" together and each data 
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point was given a time relative to time zero at the center 
of each transit event. This was done separately for the 
19-day and 39-day signals. We have been careful not to 
use a full cycle of out-of-transit data to avoid using any 
photometric measurements more than once in the input 
light curve. 

For the 1.6-day signal that repeats at regular intervals 
(since it shows no TTVs), we created a light curve by 
simply masking out the transits at the other two periods. 

3.2. Additional observations for false positive rejection 

The photometric aperture of Kepler is typically a few 
pixels across, with a scale of 3'.'98 per pixel (see below). 
High- resolution imaging of Kepler-9 was performed in or- 
der to identify neighboring stars that might be eclipsing 
binaries blended with and contaminating the target pho- 
tometry. Images were recorded with the guider camera 
of the High Res olution Echelle Spectrometer (HIRES; 
IVogt et ailll994l) on the Keck I telescope on Mauna Kea, 
in unfiltered light. The nominal sensitivity of the CCD 
from 400 to 800 nm yields an effective passband similar to 
the Kepler passband. The field of view was 43" x 57", and 
the pixel scale was 0"30 per pixel. One of these frames 
appears in Figure [1] and shows at least four stars in the 
field of view within 15" of the target. Some of these stars 
are listed in various astrometric and photometric cata- 
logs. The brightness of these companions relative to the 
target was measured using aperture photometry on four 
separate Keck images, and ranges from Am = 2.6 to 5.9 
mag. 
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Fig. 1. — Image of Kepler-9 from the HIRES guider camera on 
the Keck I telescope, obtained in seeing of 0'.'9 and clear skies. 
Companions within 15" are labeled as in Table [T] The scale of the 
image is 0''30 pix - 1 . Also indicated are the optimal photometric 
aperture (darker gray area of 8 pixels, used to extract the Kepler 
photometry) and the target aperture mask (lighter gray area of 31 
pixels, used to measure centroids) for Kepler quarter 3. 
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Fig. 2. — Sensitivity to faint companions near Kepler-9 from our 
imaging observations. Any companions above the curves are bright 
enough to be detected. J and K s limits are from AO observations 
at the Palomar 200-inch telescope, and R is from speckle observa- 
tions using the WIYN 3.5 m telescope. Companions to the right 
of the vertical dotted line at 0'.'74 cannot be responsible for the 
1.6-day signal, as they would have induced centroid motion that is 
not observed. Stars detected in our imaging observations (Table [TJ 
are marked with asterisks at their measured angular separations 
and magnitude differences in the Kepler passband. 

Speckle observations of Kepler-9 were carried out on 
2010 June 18 with the WIYN 3.5 m telescope located on 
Kitt Peak. They were taken with a two-color EMCCD 
speckle camera using narrow-band filters 40 nm wide cen- 
tered at 562 nm and 692 nm. We refer loosely to these 
passbands as V and R. The native seeing was 0"7. No 
companions with Am < 3.25 mag (R band) are present 
in the field of view centered on the target out to 1'.'8, at 
the 5-ct confidence level. Inside of 0"2 the sensitivity is 
reduced, but still allows to rule out brighter companions 
down to the diffraction limit of 0"04-0"05 (see Figure EJ. 
Details of the follow-up speckle observations in the con- 
te xt of the Kepler Mi ssion are described in more detail 
bv lHowell etaLl d2010h . 

Additionally, Kepler-9 was observed on 2010 July 
2 at the Palomar Hale 200-inch telescope with the 
near-infrared adaptiv e optics (AO) PHARO instrument 
(|Havward et alJl200l . a 1024 x 1024 Rockwell HAWAII 
HgCdTe array detector. Observations were made in the 
J (1.25 /im) and K s (2.145 /im) bands. The field of view 
was approximately 20" x 20", and the scale was 25.1 mas 
per pixel. The AO system guided on the primary target 
itself, and produced Strehl ratios of 0.05 at J and 0.3 at 
K s . The central cores of the resulting point spread func- 
tions had widths of FWHM = 0"12 at J and FWHM 
= CK'10 at K s . The closer of the companions seen ear- 
lier in the Keck images were easily detected, and we list 
them all in Table Q] along with relative positions (an- 
gular separations and position angles), relative bright- 
ness estimates, and other identifications. The sensitivity 
to faint companions was studied by injecting artificial 
stars into the image at various separations and with a 
range of Am, and then attempting to detect them both 
by eye and with an automated IDL procedure based on 
DAOPHOT. For firm detection we required the artificial 
stars to be present in both passbands. The sensitivity 
curves as a function of angular separation are shown in 
Figure [2l along with the i?-band sensitivity estimated 
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TABLE 1 

Companions to Kepler-9 identified in our imaging observations. 



Identification 


SDSS coordinates 


P 


P.A. 


AJ 


AK 3 


AKp 




(J2000) 


(") 


(deg) 


(mag) 


(mag) 


(mag) 


Kepler-9 a 


19:02:17.76 +38:24:03.2 












Comp l b 


19:02:17.91 +38:24:05.4 


2.85 


37.9 


6.84 


6.84 


5.5 


Comp 2 


19:02:18.27 +38:24:02.8 


6.04 


91.7 


4.52 


4.17 


3.7 


Comp 3 


19:02:17.29 +38:23:57.1 


8.03 


221.8 


6.25 


6.04 


5.9 


Comp 4 C 


19:02:17.69 +38:24:13.4 


10.21 


355.6 


3.59 


3.01 


2.6 



a Target is also known as 2MASS 19021775+3824032 and KIC 3323887. 

b This companion is not listed in the SDSS catalog; the coordinates are inferred from 

its position relative to Kepler-9. 

c Also known as 2MASS 19021769+3824132 and KIC 3323885. 



from the speckle observations. 

Much fainter stars with Am > 9 near a Kepler target 
could in principle be detected by examining images from 
the Palomar Observatory Sky Survey, which date back 
more than 50 years, provided the proper motion of the 
target is large enough to have shifted it by several arc 
seconds over that period. This is not possible for Kepler- 
9, since it s total proper motion a s reported in the UCAC2 
Catalog (jZacharias et a l. 2004) is only 13.7 mas yr . 
The likelihood of such faint close-in companions must 
therefore be addressed statistically, if need be. 

While the AO and speckle observations rule out the 
presence of bright neighboring stars as close as O'.'l or 
slightly less, further limits on even tighter companions 
can be placed by the spectroscopic observations ob- 
tai ned with HIRES on the Keck I telescope, described 
bv lHolman et al.l(|2010l) . since those stars would fall well 
within the 0'.'86 slit of the spectrograph. We performed 
simulations in which we added the spectrum of a faint 
star to the original Kepler-9 spectra, over a range of rela- 
tive brightnesses, and attempted to detect these artificial 
companions by examining the cross-correlation function. 
We estimate conservatively that any such stars with rel- 
ative fluxes larger than about 10-15% (Am less than 2- 
2.5 mag) would have been seen, unless their spectral lines 
are blended with those of the target. The sharp lines of 
Kepler-9, with a measured rotational broadening of only 
vsm i — 1.9 ± 0.5 kms -1 , make this rather unlikely. 

3.3. Centroid analysis 

Thanks to the very high astrometric precision of Ke- 
pler, an analysis of the motion of the photocenter of a 
target provides an effective way of identifying false pos- 
itives that are caused by background eclipsing binaries 
falling wit hin the aperture. The principles have been ex- 
plain e d by iBatalha et all (|2010|) (see also iJenkins et al.l 
l2010t iMonet et al.1 120101 ). The centroid measurements 
described below use data from quarter 3 only. In quar- 
ter 1 the Kepler-9 aperture was determined to be too 
small to optimally capture its flux, and was subsequently 
enlarged. In quarter 2 Kepler experienced undesirable 
pointing drift, which was later resolved. These problems 
complicate the centroid analysis for quarters 1 and 2, al- 
though the results are broadly consistent with the more 
reliable ones from quarter 3 presented here. 

We describe first the use of difference image analy- 
sis to demonstrate that the transit sources for all three 
Kepler-9 planets and candidates are restricted to being 
very near the target star. A difference image is formed by 
averaging several exposures near, but outside of a tran- 



sit and subtracting from this the average of all available 
exposures near transit center. This results in a typically 
isolated signal, a positive intensity with the shape of the 
point spread function (PSF) at the true spatial location 
of the transit source, and an amplitude equal to the pho- 
tometric transit depth times the direct image intensity 
for the target. Adopting 40 independent transits of KOI- 
377.03 in quarter 3 (avoiding those shortly after major 
disturbances such as a safing event, and avoiding any that 
overlap with 'b' and 'c' transits), each formed with six 
symmetrically placed exposures outside of transits (after 
a two exposure gap) and three near transit minimum, 
results in a 14-cr signal in the difference image. The cor- 
responding direct image is formed as the average of both 
in- and out-of-transit sets such that the direct and dif- 
ference images are sums and differences of precisely the 
same exposure sets. 

For Kepler-9 b four transits were used from quarter 3 
with five exposures in transit, a gap of three, then five 
more exposures on each side for out-of-transit. Kepler-9 c 
used two transits with seven exposures in-transit, a gap 
of three, and seven symmetrically placed out-of-transit 
exposure blocks. By using only exposures pulled close 
in time, and symmetrically with respect to the transits 
in use to form a difference image, this effectively im- 
poses a de-trending and avoids any complications from 
drifts on time scales longer than the average spread of 
the out-of-transit sets, which for Kepler-9 c (the widest) 
is about 9 hours. Inspection of the difference images in 
Figure [3] shows that the transit sources f or the confirmed 
'b' and 'c' planets (|Holman et al.l [20101 ) and the candi- 
date KOI-377.03 must arise from close to the target star, 
with offsets approaching one pixel easily ruled out by in- 
spection. A weighted PSF fit (or more properly, a P ixel 
Response Function (PRF) fit; see lBrvson et alJlSoiOD to 
each of the direct and three difference images of Figure [3] 
is formed using only the central 3x3 pixel area centered 
on the brightest pixel. This leads to offsets with respect 
to 'b', 'c', and KOI-377.03 of 0.007, 0.035, and 0.047 pix- 
els, respectively. For KOI-377.03 the formal error from 
a weighted least squares fit is 0.062 pixels. We have fur- 
ther assessed the errors by generating a large number of 
independent realizations of a transit signal of the KOI- 
377.03 relative intensity centered on the target coordi- 
nates. This leads to an rms scatter of 0.062 pixels. The 
noise had been increased by a factor of 1.2 in the differ- 
ence image beyond direct Poisson plus readout noise es- 
timates in order to yield this congruence of least-squares 
errors and scatter in simulations. The distribution of off- 
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Direct image Difference image Difference image Difference image 
Kepler-9 b Kepler-9 c KOI-377.03 




Fig. 3. — Direct and difference images for Kepler-9. The four panels from left to right show 128 (row) by 30 (column) regions corresponding 
to the direct and difference images for planets 'b' and 'c', and candidate KOI-377.03. The pixels returned for all stars in this area have been 
mapped into original row and column locations on the detector. Over 90% of the image area is unfilled since Kepler returns only postage 
stamps on stars of interest. The target (KIC 3323887) is indicated with circles in each panel. The locally brightest pixel is always at column 
1100 and row 273, and each display panel has been normalized by the sum of counts within the 3x3 pixels centered on [1100, 273]. The 
display range is —0.03 to 0.3. The difference images were created to isolate the signals for transits 'b', 'c', and 377.03 respectively. Most 
stars, not having variations synced with these phases, effectively disappear in the difference images. For each of the three sets of transits 
the difference image in the 3x3 pixel core appears nearly identical to the direct image, demonstrating that the true transit source must be 
near the target to a small fraction of a pixel. The difference images also reflect the expected count levels for the source to be coincident 
with the target. 



sets follows expected Gaussian statistics, e.g., in the 7472 
trials for KOI-377.03 the extreme offset is 4.2c compared 
to the expected 4.0. We have also shown that simulating 
transit signals at 0.5 and 1.0 pixel offsets from the target 
results in similar and smaller statistical scatter, respec- 
tively, as less Poisson noise is under the transit image. 
We take the scatter of 0.062 pixels to generate a 3-cr error 
circle of 0.186 pixels, or 0'.'74. This is the minimum ra- 
dius within which background eclipsing binaries cannot 
be safely ruled out from centroid analysis of the Kepler 
data itself. To place this in perspective: centroid anal- 
ysis has ruled out 98.6% of the area within the 8-pixel 
optimal aperture (>99.6% of the 31-pixel mask) of Fig- 
ure Q] as the location of potential background eclipsing 
binaries creating the KOI-377.03 signal. 

The quantitative results for all three transit sets are 
given in Table [2] Kepler-9 b shows an offset of 0.0074 
pixels between the difference and direct image relative to 
a 1-cr error of 0.0049 pixels. A 3-cr error circle in which 
background binaries cannot be excluded from the cen- 
troid analysis of Kepler data itself is only 0'.'06. Kepler- 
9 c is the only case of the three showing a formal incon- 



sistency with the offset being 5cr; however, even if we 
combine the offset and 3-cr formal error any background 
eclipsing binaries outside of a radius of 0'.'22 are excluded 
as the transit signal source. Clearly for all three tran- 
sit sets, with the 3-cr error circles comfortably under 1" 
all of the known companions from high-resolution imag- 
ing shown in Table Q] are safely ruled out as sources of 
the photometric transit signal. It is worth noting that 
the formal (and equal to scatter from Monte-Carlo sim- 
ulations) error on radial offsets is approximately equal 
in pixel units to the inverse photo metric sign al-to-noise 
(S/N) ratio, as expected (see, e.g.. lKinglll983l ). 

Further confirmation that at least the two deeper sig- 
nals seen in Kepler-9 are not due to known stars in the 
scene can be obtained by placing simulated eclipses on 
the known stars in the aperture, and comparing them 
with the observations. The scene in the aperture is 
modeled using star s in the Kepler Input Catalog (KIC; 
Lath anTet al.l l2005l ). supplemented by the stars in Ta- 
ble [1] All stars within a PRF size (15 pixels) in row or 
column of Kepler-9 's aperture are included. To generate 
the modeled out-of-transit image, the measured PRF is 
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TABLE 2 

PRF CENTROID MEASUREMENTS ON KEPLER-9 DIRECT AND 
DIFFERENCE IMAGES. 



Type 


fntcnsity 


Column 


Row 


Radius 




(«-) 


(pix) 


(pix) 


offset (pix) 


Kcpler-9 b 


Direct 


5.170 x 10 7 


1099.6989 


273.4557 




Difference 


4.360 x 10 5 


1099.7058 


273.4584 




S/N ; Offset 


174 


0.0069 


0.0027 


0.0074 


Errors 


2506 


0.0029 


0.0039 


0.0049 


Kepler-9 c 


Direct 


5.132 x 10 7 


1099.6999 


273.4654 




Difference 


3.833 x 10 5 


1099.7107 


273.4990 




S/N ; Offset 


128 


0.0108 


0.0336 


0.0353 


Errors 


2991 


0.0039 


0.0056 


0.0068 


KOt-377.03 


Direct 


5.401 x 10 7 


1099.6946 


273.4408 




Difference 


1.493 x 10 4 


1099.7131 


273.4836 




S/N ; Offset 


14.3 


0.0185 


0.0428 


0.0047 


Errors 


1046 


0.0351 


0.0511 


0.0062 



Note. — The first two lines of each block present Intensity and 
two coordinate position PRF fit results for the Direct and Difference 
images, respectively. The third line shows the photometric signal-to- 
noisc for the intensity in the Difference image, then the offset in posi- 
tion of the preceding two lines, with the last entry being the quadra- 
ture sum of the column and row offsets. Errors refer to the PRF fit 
to the Difference image. The scale is 1 pixel — 3. 98. 

placed at each star's location on the focal plane, scaled 
by that star's flux. This provides the contribution of each 
star to the flux in the aperture's pixels. For each star s$ 
in the aperture, the depth d Si of a transit is computed 
that reproduces the observed depth in the aperture pix- 
els. An in-transit image for each Sj is created as in the 
out-of-transit image, but with the flux of Sj suppressed 
by 1 — d Si . These model imag es are subject to errors in 
the PRF (|Brvson et al.l I201QB . so they will not exactly 
match the sky. 

A flux-wcighted centroid is computed for the out-of- 
transit image and the in-transit-image generated for each 
star in the aperture. This produces row and column cen- 
troid offsets AR and AC, and the centroid offset distance 
D = VAR 2 + AC 2 . 

To compare these modeling results with observation we 
must make low-noise centroid measurements from the ob- 
served pixel data. We do this by creating out-of-transit 
and in-transit images from de-trended, folded pixel time 
series. For each pixel time series, the de-trending op- 
eration has three steps: 1) removal of a median-filtered 
time series with a window size equal to the larger of 48 
cadences or three times the transit duration; 2) removal 
of a robust low-order polynomial fit; and 3) the appli- 
cation of a Savitzky-Golay filtered time series of order 
3 with a width of 10 cadences. The Savitzky-Golay fil- 
ter is not applied within 2 cadences of a transit event, 
so the transits are preserved. The resulting pixel time 
series are folded with the transit period. Each pixel in 
the out-of-transit image is the average of 30 points taken 
from the folded time series outside the transit, 15 points 
on either side of the transit event. Each pixel in the 
in-transit image is the average of as many points in the 
transit as possible: seven for Kepler-9 b and Kepler-9 c, 
and four for KOI-377.03. Centroids are computed for the 
in-transit and out-of-transit images in the same way as 



the modeled images. 

Uncertainties of these centroids are estimated via 
Monte Carlo simulation, where a noise realization is in- 
jected into 48-cadence smoothed versions of the pixel 
time series for each trial. A total of 2000 trials are per- 
formed each for Kepler-9 b, Kepler-9 c, and KOI-377.03. 
The in- and out-of-transit images are formed using the 
same de-trending, folding and averaging as the flight 
data. The measured uncertainties are in the range of 
a few times 10~ 5 pixels. 

Table [3] shows the resulting measurements of the cen- 
troids from quarter 3 pixel data, along with the Monte- 
Carlo-based 1<T uncertainties. The centroids are con- 
verted into centroid offsets and offset distance with prop- 
agated uncertainties. Table [4] shows the offset distance 
D predicted by the modeling method described above 
for each target in the aperture. We see that when the 
transit is on Kepler-9 itself we expect a measurable cen- 
troid shift for Kepler-9 b and Kepler-9 c. In this case 
the modeled centroid shift is about 3.7 times larger than 
that observed, though the signs of the offsets agree. This 
exaggeration of the centroid offset has been traced to in- 
accuracies in the KIC used to create the model images. 
Therefore, the centroid shifts in Table 0] should be scaled 
by a factor 1/3.7. If the transit were on one of the com- 
panion stars in the aperture, then the modeled centroid 
shift would be an order of magnitude larger than ob- 
served for Kepler-9 b and Kepler-9 c, ruling out the com- 
panion stars as the source for these signals. Companion 
stars are not as definitively ruled out for the KOI-377.03 
transit by this technique. After scaling the centroid off- 
sets as above, modeled transits on companions 3 and 4 
have offsets that are are about 2.5cr, while companion 2 
is 1.9<7 and companion 1 is less than ler. The modeled 
transit on Kepler-9, however, is much smaller, consistent 
with the observed transit offset for KOI-377.03. 

3.4. BLENDER analysis of Kepler-9 b and c 

As an initial test, we modeled the light curves for each 
of these two signals assuming they are the result of an 
eclipsing binary physically associated with the target, 
i.e., at the same distance (hierarchical triple). For this 
case the isochrone for the binary was taken to be the 
same as that of the primary, and corresponds to [Fe/H] 
= +0.12 and an age of 1 Gyr. The secondary and tertiary 
masses were allowed to vary freely between 0.10 Mq (the 
lower limit in the models; see footnote |22|) and 1.40 Mq, 
as mentioned earlier, seeking the best fit to the photom- 
etry. The inclination angle was also free, and the orbits 
were assumed to be circular. In both Kepler-9 b and c, 
which have similar transit signals, we find that the best 
fitting hierarchical triple blend model corresponds to sec- 
ondaries that are approximately 1.0 and 0.5 mag fainter 
than the primary respectively, and tertiaries that are at 
the lower limit of the isochrone range (late M dwarfs). 
However these fits give a poor match to the photometry: 
BLENDER is unable to simultaneously reproduce the total 
duration of the transit and the central depth, given the 
constraints on the brightness and size of the stars from 
the isochrones. This type of blend scenario is therefore 
clearly ruled out. We illustrate this for Kepler-9 b in Fig- 
ure 01 where the best-fit planet model is also shown for 
reference. Much better matches to the data can be found 
if additional light from a fourth star along the line of sight 
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TABLE 3 

Observed centroid shifts for Kepler-9 b, Kepler-9c, and KOI-377.03. 





Kepler-9 b 




Kepler-9 c 




KOI-377.03 




AR 


2.52 x 10~ 4 ± 7.78 x 10 


-5 


1.65 x 10~ 4 ± 9.34 x 10 


-5 


-1.24 x 10" 7 ± 6.19 x 10" 


■ r > 


AC 


-2.23 x lO" 4 ± 7.55 x 10 


-5 


-2.40 x 10" 4 ± 9.21 x 10 


-5 


8.39 x 10" 6 ± 5.82 x 10" 


5 


D 


3.41 x 10~ 4 ± 7.66 x 10 


-5 


2.91 x 10~ 4 ± 9.25 x 10 


-5 


8.39 x 10" 6 ± 5.82 x 1(T 


5 


D/a 


4.44 




3.15 




0.14 




Note. 


— The measurements are g 


iven 


in pixel units, and the scale 


is 3'. 


'98 per pixel. 





TABLE 4 

Modeled centroid shifts due to transits on the known stars in the aperture with depths 

that reproduce the observed depth. 



Object 


Modeled D 


D/a 


Object 


Modeled D 


D/a 


Object 


Modeled D 


D/a 


Kepler-9 b 


1.22 X 10" 3 


16.0 


Kepler-9 c 


1.10 X 10" 3 


11.9 


KOI-377.03 


5.02 x 10~ 5 


0.86 


Comp 1 


depth > 1 




Comp 1 


depth > 1 




Comp 1 


1.45 x 10~ 4 


2.49 


Comp 2 


9.77 x 10" 3 


127 


Comp 2 


8.75 x 10~ 3 


94.6 


Comp 2 


4.01 x 10~ 4 


6.89 


Comp 3 


depth > 1 




Comp 3 


depth > 1 




Comp 3 


5.33 x 10~ 4 


9.16 


Comp 4 


1.44 X 10" 2 


188 


Comp 4 


1.29 X 10~ 2 


139 


Comp 4 


5.90 x 10~ 4 


10.1 


Note. — 


Shifts are given 


in pixel units, and 


the scale is 3. 


38 per pixel. For Kepi 


er-9 b and Kcp] 


er-9 c 



transits on some companions can be ruled out because they require depth > 1. 



is incorporated into the model, providing extra dilution. 
We find that this fourth star is required to be nearly as 
bright as the primary, and the optimal model changes in 
such a way that the secondary also becomes as bright as 
the primary (so that its size enables the duration of the 
transits to be reproduced), while the tertiary remains a 
small star. This rather contrived scenario requiring two 
bright stars that are nearly identical to the main star 
would be easily recognized in our high-resolution imaging 
for separations larger than about O'.'l (see, e.g., Figure[5]), 
in our centroid analysis for separations larger than 0'.'06, 
or would otherwise produce obvious spectroscopic signa- 
tures unless all three bright objects happened to have the 
same radial velocity. 

We next considered blends involving eclipsing bina- 
ries in the background, by removing the constraint on 
the distance. In this case a solar-metallicity isochronc 
was adopted for the binary, with a representative age for 
the field of 3 Gyr. We explored a wide range of relative 
distances, and we first considered main-sequence stars 
only, again with circular binary orbits. The results for 
Kepler-9 b and c are once again similar to each other, 
and we illustrate them for Kepler-9 c in Figure [5] The 
axes correspond to the distance modulus difference AS 
as a function of the tertiary mass. Contours represent 




Planet model 
Blend model 



-0.010 -0.005 0.000 0.005 0.010 

Orbital phase 

Fig. 4. — Light curve of Kepler-9 b (P = 19.24 days) with the 
best fit blend model for the case of a hierarchical triple (candidate 
+ physically associated eclipsing binary) . The best fit planet model 
is shown for reference. The poor fit of the blend model rules out 
this configuration. 



constant differences in the % 2 of the fit compared to the 
best-fit planet model, and are labeled in units of the sta- 
tistical significance of the difference (a). We draw two 
main results from this figure. One is that the light curve 
fits strongly prefer the smallest available tertiary masses 
from the isochrones (0.10 Mq), and would in fact yield 
better fits for even smaller tertiaries (i.e., planets). Addi- 
tionally, the best solutions cluster toward equal distances 
for the binary and the primary star, effectively converg- 
ing toward the equivalent of the hierarchical triple sce- 
nario considered earlier. No acceptable solutions exist 
with the binary at a significant distance behind the pri- 
mary star. The best fit to the light curve of Kepler-9 c 
is similar to the one shown in Figure [4] (dashed curve) , 
which is not particularly good. The AS vs. tertiary mass 
diagram for Kepler-9 b is qualitatively the same. Allow- 
ing the secondary to be a giant star gives a very poor 
fit to the photometry: the duration of the transit is very 
much longer than observed, there is out-of-eclipse mod- 
ulation due to distortions in the giant, and all solutions 
place the binary at an implausibly large distance. We 
conclude that blend configurations involving background 
eclipsing binaries in which the tertiary is a star are not 
a viable explanation for either of these two signals. 

We then explored background eclipsing binaries in 
which the tertiaries are planets rather than stars. This 
allows their radii to be smaller, possibly providing a bet- 
ter fit to the Kepler photometry. The orbits were consid- 
ered to be circular, as before. Figure [5] shows the results 
for Kepler-9 c, this time in the plane of separation versus 
secondary mass. Once again the fits tends to favor an 
equal distance for the binary and the primary star, and 
background scenarios with the binary far behind provide 
unacceptably poor matches to the light curve. A second 
noteworthy result is that these solutions have a strong 
preference for secondary stars that are quite similar to 
the primary. All acceptable fits to the light curve corre- 
spond to relatively bright secondaries with AKp < 1.5 
mag (see Figure [6|). The best of these solutions is of 
about the same quality as a planet model, and has a sec- 
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0.0 0.2 0.4 0.6 0.8 1.0 

Tertiary mass (M 6 ) 
Fig. 5. — Map of the x 2 surface (goodness of fit) corresponding 
to a grid of blend models for Kcpler-9c (P = 38.91 days) involving 
background eclipsing binaries with circular orbits. The separation 
between the binary and the primary is expressed in terms of the 
distance modulus difference. Contours are labeled with the \ 2 
difference from the best planet model fit (expressed in units of 
the significance level of the difference, <r), and are plotted here as 
a function of the mass of the tertiary star. The dashed line at 
O.lMg indicates the lower limit to the tertiary mass set by the 
model isochrones we use. 

ondary of mass 0.98 Mq that is only 100 K cooler and 
0.3 mag fainter than the primary in the Kepler band. 
This somewhat artificial case of "twin" stars is a result 
we have seen often in simulations for other Kepler candi- 
dates. The tertiary in this type of blend solution comes 
out about V2 larger than in a planet model because the 
transit is diluted by another star of approximately equal 
size and brightness. One may debate whether this situa- 
tion should actually be referred to as a "false positive" for 
Kepler-9, since the signal would still correspond to a gas 
giant planet, only that this planet would be V2 larger, 
and it would be orbiting a different star. Alternatively, 
it could be thought of simply as an overlooked dilution 
factor in a true planetary system. In any event, the lack 
of evidence for this bright twin star in the spectroscopy 
or in our high-resolution imaging or centroid analysis for 
Kepler-9 does not support this scenario. 

As a particular case of this family of configurations, 
we examined blends in which the star-planet pair is con- 
strained to be at the same distance as the primary, i.e., 
effectively in a hierarchical system. The secondary prop- 
erties were therefore taken from the same isochrone as 
the primary, and the orbits were assumed to be circular. 
An excellent fit to the light curve is possible for a tertiary 
that is about \[2 larger than in a true planet model, but 
not surprisingly, we find once more that the secondary 
must be as bright as the primary. 

Additional tests were run to examine the impact of 
changing the age adopted for the isochrone of the sec- 
ondary in a background star-planet pair, or the addition 
of light from a fourth object in the aperture. In the first 
case, changing the age from 3 Gyr to 1 Gyr produced a 
small shift of the contours in Figure |5] downward and to 
the right that is simply due to the change in intrinsic 
brightness of the secondary star, and does not alter our 
conclusions. Adding "fourth light" further attenuates the 
eclipses of the star-planet pair. To compensate, BLENDER 




0.2 0.4 0.6 0.8 1.0 1.2 1.4 
Secondary mass (M ) 
Fig. 6. — Map of the x 2 surface corresponding to a grid of blend 
models for Kepler-9 c involving background eclipsing systems in 
which the tertiary is a (dark) planet, in a circular orbit around the 
secondary. Contours are labeled with the \ 2 difference from the 
best planet model fit (expressed in units of the significance level of 
the difference, <r). Two dashed lines of equal magnitude difference 
(AKp) are indicated, and show that all viable blend fits (with 
confidence level < 3a") have secondaries that are bright enough to 
have been detected spectroscopically (AKp < 2). 

requires a slightly deeper eclipse, and in order to pre- 
serve the shape of the signal (total duration, and slope 
of ingress/egress), this is achieved by bringing the sec- 
ondary closer to the primary. As a result, for relatively 
bright fourth light the contours are shifted downward by 
approximately the difference in magnitude between the 
primary and the fourth star, again without changing the 
conclusions. 

One may also imagine blend scenarios in which the 
eclipsing binary is in the foreground, rather than the 
background. We explored this possibility by extending 
the simulations to negative values of Ai5. As before, we 
adopted circular orbits and a 3 Gyr isochrone for the fore- 
ground system. Binaries with stellar tertiaries are clearly 
ruled out as they yield fits to the light curve that do not 
match its shape, and additionally they predict a fairly 
obvious secondary eclipse that is not seen in the data. 
We focus therefore on blends in which the tertiary is a 
planet, and we illustrate the results for Kepler-9 c. In 
this case we find there are many acceptable solutions 
with x 2 values differing from the best planet fit at the 
level of la or less. These solutions span a range of sec- 
ondary masses and a range of foreground separations, 
implying a wide range not only in apparent brightness 
for the secondary, but also in color. Models in which the 
secondaries are brighter than the primary and of signifi- 
cantly different spectral type would be inconsistent with 
the spectroscopic parameters derived for Kepler-9, and 
are excluded. Plausible solutions remain, in principle, 
for fainter foreground secondaries, which necessarily in- 
volve later-type stars. We find that of these, the only 
ones that yield acceptably good fits to the Kepler pho- 
tometry, with \ 2 values differing from the planet fit by 
less than 3rj, correspond to secondaries that are within 
about 1.5 mag of the the primary in brightness, and are 
of course redder. These would be valid blend configu- 
rations so long as the secondaries are close enough to 
the primary to be spatially unresolved (angular separa- 
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tions < OVl), and at the same time faint enough to have 
gone undetected in the spectra. Stars that are within 
~2-2.5 mag of the primary would generally have been 
seen spectroscopically, as indicated in Sect. 13.21 and this 
would exclude these remaining foreground blend config- 
urations. Nevertheless, to be conservative, let us assume 
for the moment that a star 1.5 mag fainter than the pri- 
mary has still managed to elude detection in our spec- 
tra. This corresponds to the faintest secondary in a fore- 
ground blend scenario that still allows for a satisfactory 
fit to the light curve, and would be the most difficult case 
of this kind to disprove. This fit is shown in the top panel 
of Figure and is statistically indistinguishable from a 
planet model fit. The secondary in this configuration is 
an M2 dwarf (M = 0.56 Mq) 1.53 mag fainter than the 
primary, eclipsed by a 0.91i?j up planetary companion, 
and is located at a distance of 300 pc. The primary in 
this scenario is at 750 pc. 

Other properties of this particular blend such as mag- 
nitudes and colors can be computed easily with BLENDER, 
and compared with observations. Apparent magnitudes 
for Kepler-9 are available from the KIC for a variety 
of passbands including Sloan griz, a special-purpose 
passband referred to as D51 (centered on the Mg lb 
triplet at 518.7 nm), and JHK S from the 2MASS cat- 
alog. The lower panel of Figure [7] shows various color 
indices (Kp — A) predicted by BLENDER both for the pri- 
mary star alone and for the blend. Those of the primary 
are well reproduced by the model, and we find that a 
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Fig. 7. — Top: Light curve of Kepler-9 c with the best fit blend 
model for the case of contamination by a foreground eclipsing pair 
with a circular orbit in which the tertiary is a planet. The pair 
consists of an M2 dwarf (0.56 Mq, 0.58 Rq) and a 0.91 -Rj up com- 
panion 450 pc in front of the primary, which is at 750 pc. This 
fit is statistically indistinguishable from best fit planet model, also 
shown for reference. Bottom: Measured colors for Kepler-9 (dots) 
compared with the predictions from the blend model in the top 
panel. A small amount of extinction (0.15 mag kpc -1 ) has been 
included in these predictions. The results without considering ex- 
tinction differ little, and are shown with dotted lines. The color 
measurements clearly rule out such a blend. 



small amount of interstellar extinction leads to an even 
better match (solid line in the figure). The colors of the 
blend, on the other hand, disagree with the measured 
colors, and deviate by more than 0.4 mag for the reddest 
index, Kp — K s . We are therefore able to exclude, solely 
on the basis of its color, this most difficult of the sce- 
narios involving foreground star-planet pairs that could 
mimic the 19-day and 39-day signals in the light curve of 
Kepler-9. Larger-mass secondaries would not be as red 
and still allow for good fits to the photometry, but they 
are intrinsically brighter and would be recognized more 
easily. 

The above simulations have all assumed circular orbits 
for the blended eclipsing binaries or star-planet pairs, 
which is not necessarily realistic given the relatively long 
periods of Kepler-9 b and c. Eccentricity affects the speed 
of the secondary and tertiary in their relative orbit, and 
therefore can change the duration of the transit, mak- 
ing it shorter or longer than in a circular orbit, depend- 
ing on the orientation (longitude of periastron, uS). It 
also changes the impact parameter, all else being equal. 
And finally, it shifts the location of the secondary eclipse. 
The magnitude of these effects is illustrated in Figure [8] 
for eccentricities between 0.1 and 0.7. The most im- 
portant effect for our purposes is on the transit dura- 
tion. Given a fixed (measured) duration, eccentric or- 
bits may allow blends with smaller or larger secondary 
stars than in the circular case to still provide satisfactory 
fits to the light curve, effectively increasing the pool of 
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Fig. 8. — Effect of eccentricity on the duration of transits rela- 
tive to the circular orbit case (A/A c ; rc ), on the impact parameter 
(b/6 c i rc ), and on the displacement (</> scc — 0.5) of the secondary 
eclipses relative to phase 0.5, all shown as a function of the longi- 
tude of periastron u>. The different curves correspond to eccentric- 
ities from 0.1 to 0.7 in steps of 0.1. 
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Fig. 9. — Same as Figure[6](Kepler-9c), restricted to star-planet 
orbits having e = 0.3 (concentration of contours on the left) and 0.5 
(right), and uj = 90°. This orientation corresponds to transits that 
occur at periastron. Comparison with Figure [6] shows that these 
solutions allow for more massive (larger) secondary stars than in 
the case of circular orbits, but the brightness of these blends is still 
within 2 magnitudes of the target, and is ruled out by spectroscopy. 



potential false positives. The limiting cases correspond 
to to = 90° and 270°, in which the line of apsides is 
aligned with the line of sight and the transit occurs at 
periastron (accommodating larger secondaries) or apas- 
tron (smaller secondaries), respectively. Extensive simu- 
lations for these two extreme situations show that allow- 
ing for eccentric orbits does not change our conclusions 
regarding hierarchical triple systems, background eclips- 
ing binaries, or background star-planet scenarios. We 
show this for the latter blend category in Figure [9j illus- 
trated for the case of orbits with eccentricities of 0.3 and 
0.5, and uj = 90°. Comparison with Figure |6] indicates 
that in both cases the blends are still bright enough that 
we would have seen signatures of them in the spectra of 
Kepler-9. Larger eccentricities of e = 0.7 result in secon- 
daries that are brighter still. For eccentric orbits oriented 
such that transits take place at apastron (uj = 270°), we 
only find acceptable fits to the light curves for eclipsing 
star-planet pairs that are in the foreground (and involve 
smaller stars). However, as was the case for circular or- 
bits, those blends are either too bright, too red, or both, 
and are thus also excluded. 

The above, fairly exhaustive exploration of parame- 
ter space with BLENDER allows us to conclude that no 
configuration involving an eclipsing binary (or an eclips- 
ing star-planet pair), either in the foreground or in the 
background, is able to provide a reasonable explanation 
for the signals of Kepler-9 b and c (see Table [5] for a 
summary of the configurations tested, and the results). 
Many scenarios lead to light curves that match the de- 
tailed shape of the transit events, but none are simul- 
taneously consistent with all of the other observational 
constraints. This includes spectroscopy, high- resolution 
imaging, centroid measurements, and photometry (col- 
ors). Therefore, even ignoring the evidence from TTVs, 
these results fully support the planetary nature of these 
objects and demonstrate the usefulness of BLENDER for 
validating transiting planet candidates from Kepler. 



3.5. BLENDER analysis of KOI-377. 03 

We proceed next to examine false positive scenarios for 
the shallowest signal in Kepler-9, with P — 1.59 days, 
which would correspond to a super-Earth-size planet. 
Because the period is so short in this case, and tidal 
forces in such bin ary systems have likely circularized 
the orbit (see, e.g., Mazch 2008, and references therein), 
we do not consider non-zero eccentricities. Additionally, 
blends in which the secondary star is a giant need not be 
considered, as those cases are obviously ruled out because 
of the short orbital period and small implied semimajor 
axis of the orbit. 

As for the larger signals considered above, hierarchi- 
cal triple systems in which the tertiary is a star fail to 
provide good fits to KOI-377. 03. A good match to the 
Kepler photometry can be found when the tertiary is al- 
lowed to be a much smaller object (i.e., a planet), but as 
was the case earlier, it requires a secondary that is very 
similar to the primary in brightness. The resulting size of 
the eclipsing object is larger than in a planet model, 
or slightly over 2 . This type of configuration was 
ruled out earlier based on the high-resolution imaging 
and the spectroscopy. Small tertiaries with appreciable 
mass, such as white dwarfs, induce tidal distortions on 
the primary due to the short orbital period that lead to 
significant out-of-eclipse variations in the light curve (el- 
lipsoidal variability). These modulations are not seen in 
the photometry for Kepler-9, and such false positives are 
therefore also excluded. 

Blends with an eclipsing binary in the background (the 
tertiary being a star) are able to reproduce the light curve 
just as well as a planet model. In Figure [TU] we illustrate 
those results by showing the area of allowed parameter 
space in a diagram of distance modulus difference as a 
function of secondary mass. Acceptable fits with \ 2 dif- 
fering from the planet model by less than 3tr are possible 
over a wide range of relative separations (4.5 < AS < 9), 
but the secondaries are restricted to a relatively narrow 
interval in mass centered on the mass of the primary. 
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Fig. 10. — Map of the \ 2 surface for KOI-377. 03 corresponding 
to a grid of blend models involving background eclipsing binaries 
with stellar tertiaries. Contours are labeled with the x 2 difference 
from the best planet model fit (expressed in units of the signifi- 
cance level of the difference, cr). The dashed lines indicate levels of 
equal apparent magnitude difference AKp between the background 
binary and the primary star. 
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TABLE 5 

Summary of blend configurations tested for. Kepler-9 b and c. 



False positive configuration 51 



Result 



Blends ruled out 



Hierarchical triple with stellar tertiary, MS 

• Circular and eccentric orbits 

• Added fourth light 



Hierarchical triple with planetary tertiary, MS 
• Circular and eccentric orbits 



Background EB with stellar tertiary 

• Circular and eccentric orbits, MS and giants . 

Background EB with planetary tertiary, MS 

• Circular and eccentric orbits 

• 1 Gyr isochrone for secondary 

• Added fourth light 



Foreground EB with stellar tertiary, MS 
• Circular and eccentric orbits 



poor fits/sec. eel. 
twin star 



twin star 



poor fits 



twin star 
little change 
little change 



poor fits/sec. eel. 



Yes 

Yes (imaging/spec. /centr.) 
Yes (imaging/spec. /centr.) 
Yes 

Yes (imaging/spec. /centr.) 
Yes (imaging/spec. /centr.) 
Yes (imaging/spec. /centr.) 

Yes 



Foreground EB with planetary tertiary, MS 
• Circular and eccentric orbits 



too bright/too red Yes (imaging/spec. /centr. /color) 



a 3 Gyr isochrone and solar metallicity assumed for background and foreground stars, unless otherwise indicated. Ab- 
breviations: MS — main sequence secondary; imaging/spec. /centr. — high-resolution imaging, spectroscopy, and centroid 
analysis; sec. eel. — secondary eclipses predicted but not observed; EB — eclipsing binary. 



These eclipsing binaries are all very distant and faint 
(Kp « 19-22), and have no effect on the colors of the 
blend. The more distant scenarios place the binary at 
implausibly large distances of up to 42 kpc (more than 
10 kpc above the Galactic plane). The nearest config- 
uration (AS = 4.5; see Figure [TU)) has the binary at a 
distance of 5.3 kpc, and the primary at ~670 pc. The 
secondary in this model is a late G star 5.5 mag fainter 
than the primary in the Kepler passband, eclipsed by 
a late M dwarf that produces no detectable secondary 
eclipse. The predicted brightness of this binary precisely 
matches that of the closest companion identified in the 
AO images (Comp 1, Tabled]), located 2"85 NE of the 
target. However, this and all wider visual companions 
are already ruled out at more than the 3-er confidence 
level by the lack of centroid motion, which would have 
revealed any blended eclipsing binaries at angular sepa- 
rations larger than about 0'.'74 (Sect. l5T3"|) . Even without 
this constraint from astrometry, the predicted J—K s color 
of the secondary in this blend is considerably redder than 
measured for this close AO companion, which would also 
disqualify it. Eclipsing binaries that are between 5 and 
~8.5 magnitudes fainter than the main star provide ac- 
ceptably good fits to the light curve (see Figure [TU]), and 
if they were angularly closer than 0'/74 from the target 
they may not be detectable in our AO or speckle obser- 
vations, in our centroid motion analysis, nor in our spec- 
tra. They remain viable blend configurations, and would 
necessarily be at distances greater than 5 or 6 kpc. An 
example is shown in Figure [TTJ to illustrate that the fit 
is indistinguishable from a planet fit. 

In the above calculations we have ignored interstel- 
lar extinction. However, given the large distances for 
the binaries in some of these blend configurations, it is 
worth exploring the effect of dust more carefully, which 
we have done by repeating the BLENDER simulations us- 
ing a representative differential extinction coefficient of 
0.5 mag kpc -1 . The results are shown in Figure [T2l and 
indicate that the blend scenarios providing good fits to 
the Kepler photometry of KOI-377.03 are systematically 



shifted to smaller distances compared to the previous cal- 
culations. Their apparent brightness, however, changes 
relatively little, as can be seen by comparing the lines of 
equal AKp with those in Figure [TUJ Therefore, the over- 
all impact of differential extinction on the permitted area 
of parameter space in terms of observable parameters is 
not as significant as might have appeared. 

Allowing the tertiary to be a smaller object such as a 
planet opens up a different area of parameter space for 
permissible background blends (Figure I13[) . When in- 
cluding extinction as before, acceptable fits to the light 
curve are possible for AS values from zero up to about 
4.5. This upper limit corresponds to distances for the 
star-planet pair of about 4.8 kpc (with AKp rj 6, or 
apparent magnitudes of Kp ss 20), and is set by the 
maximum size of 1.8i?j up we have allowed for a planet. 
As in the configurations described before, the solutions 
constrain the secondary masses to be near that of the pri- 
mary in order to match the detailed shape and observed 
duration of the transits, with a range from about 0.9 M© 
to 1.2 Mq. Therefore, the color of the blend is not as 
useful a discriminant in this case. Secondaries that are 
less than about 2 mag fainter than the primary would 
have been seen spectroscopically. This excludes a good 
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Fig. 11.— Example of a blend model fit to KOI-377.03 involving a 
background eclipsing binary with a stellar tertiary (solid line). The 
secondary is similar in spectral type to the primary and 5.2 mag 
dimmer, and the tertiary is a late M dwarf. The eclipsing pair is 
6 kpc behind the primary. This fit is statistically indistinguishable 
from the best fit planet model, which is shown with a dashed line. 
The Kepler observations have been binned for clarity. 
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Fig. 12. — Same as Figure llOl including the effect of differen- 
tial extinction in the amount of 0.5 mag kpc -1 . The net effect of 
extinction is to compress and shift the contours toward smaller rel- 
ative distances. The dashed lines indicate levels of equal apparent 
magnitude difference AKp between the background bin ary and the 
primary star, and are the same as shown in Figure [TOl 



fraction of the space of parameters, as indicated by the 
lower dashed line in Figure 1131 Of the remaining blends 
of this kind between AKp — 2 and AKp = 6, only the 
ones with angular separations smaller than about 1" are 
allowed by the constraints from our AO imaging (see Fig- 
ure [2]), but the centroid analysis is even more restrictive 
and rules out stars outside of 0'.'74. At closer separa- 
tions, the high-resolution images rule out all blends that 
are brighter than the sensitivity limit indicated in Fig- 
ured] (i.e., those that fall above the curves). 

As expected from the fixed duration and depth of 
the transit-like signal of KOI-377.03, the size of the 
tertiary in these configurations correlates with the sec- 
ondary mass. Due to this correlation, small tertiarics 
with R < 0.3.Rj U p (roughly Neptune-size, and smaller) 
are further excluded because the eclipses they produce 
are already very shallow, and further dilution by the 
primary would make them too shallow to fit the pho- 
tometry. In order to avoid this, the secondaries in those 
blends must be relatively small late-G type stars that are 
nearby, and would therefore be bright enough (AKp < 2) 
that they would have been detected in our spectra as a 
second set of lines. Thus, BLENDER effectively places lim- 
its not only on the secondary, but also on the size of 
the tertiary (see Figure fl4|) . In particular, blends with a 
white dwarf eclipsing a background star are also ruled out 
for the same reason described above. Additionally, the 
predicted light curves for such cases with white dwarf 
tertiaries show ellipsoidal variability, which is not ob- 
served. 26 Many of the larger tertiaries correspond to gas 
giants (Saturn-size, or larger), which implies a qualita- 
tive difference in their nature compared to the alternate 

26 We note, for completeness, that the mass of a white dwarf 
would generally also be sufficient to induce tidal synchronization in 
the secondary star, resulting in line broadening that could in prin- 
ciple render it more difficult to detect in the spectrum. However, 
given the orbital period and typical secondary sizes, we estimate 
the rotational broadening to be no more than ~30kms — 1 , which 
should still allow that star to be seen spectroscopically if it were 
bright enough. In any case, white dwarfs are excluded as viable 
tertiarics for the reasons mentioned in the text. 
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Fig. 13. — Same as Figure [T2l but for the case in which the 
tertiary is a planet instead of a star. Differential extinction is 
included. Kinks in the contours are an artifact of the discreteness 
of our grid. The dashed lines indicate levels of equal apparent 
magnitude difference AKp between the background secondary and 
the primary star. The lower of these lines represents the constraint 
from the spectroscopy for Kepler-9 (see text). 
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Fig. 14. — Same as Figure fl3l shown here as a function of the 
tertiary radius. Kinks in the contours and closed inner contours 
for the 1-ct and 2-cr levels are an artifact of the discreteness of 
our grid. We indicate with a dashed line the lower limit for the 
size of the tertiaries that is set by the spectroscopic constraint on 
presence of bright stellar companions (see text). The size of the 
Earth, Neptune, and Saturn are also indicated for reference. 

model of a true Earth- or super-Earth-size planet. In this 
sense these blends may properly be considered "false pos- 
itives" , as opposed to the configurations discussed earlier 
requiring twin stars, which only change the tertiary ra- 
dius by v2- 

Finally, we examine the possibility that the true period 
of the KOI-377.03 signal is twice the nominal value. Al- 
ternating events would then correspond to the primary 
and secondary eclipses of a blended eclipsing binary (the 
tertiary being a star in this case), which may in gen- 
eral be of different depth. In KOI-377.03 there is no 
compelling evidence for a depth difference between odd- 
and even-numbered events, but this is difficult to estab- 
lish in a faint star such as this for a signal that is only 
0.2 mmag deep. The results of extensive simulations with 
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Fig. 15. — Same as Figure lTOl but for the case in which the orbital 
period is assumed to be twice the nominal value (2P = 3.185702 
days). 

BLENDER for this scenario are illustrated in Figure [15l 
The best fits correspond to blended binaries far in the 
background, and do not indicate a significant difference 
in depth between the primary and secondary eclipses. 
However, these fits provide only a poor representation of 
the Kepler light curve, and can therefore be confidently 
ruled out. This is seen in Figure 1161 The top panel 
shows the closest fit to the full light curve together with 
the data, and in the bottom panel we have binned the 
measurements to facilitate the comparison. This solution 
involves an eclipsing pair of mid-G dwarfs at 21 kpc and 
7.2 mag fainter than the primary, and the fit is visibly 
worse than that corresponding to a planet at half the 
period, which is shown with the dashed line. 

In summary, the BLENDER analysis of this section cou- 
pled with constraints from spectroscopy, high-resolution 
imaging, and centroid motion measurements rules out a 
large fraction of the false positives that could produce the 
1.59-day signal, but not all (see Table [6]). The remaining 
configurations involve main-sequence background stars 
that are similar to the primary in spectral type (late F 
to early K, or about 0.9 M Q to about 1.2 Af Q ), and are 
eclipsed either by another main-sequence smaller star or 
by a planet with R p > 0.3i?j up (Neptune-size or larger). 
These blends range in apparent brightness from Kp s» 19 
to Kp « 22 for stellar tertiaries, and from Kp « 16 to 
Kp s» 20 if the tertiaries are planets, and must be closer 
than 0774 from the target. At separations under 0774 our 
imaging observations allow us to rule out the brighter of 
these blends, and only the ones with Am below the sen- 
sitivity curves in Figure [2] would remain undetected. 

3.6. Likelihood of remaining blend scenarios for 
KOI-377.03 

In the previous section we have considered a wide va- 
riety of possible blend scenarios for KOI-377.03 involv- 
ing secondaries of different spectral types (main-sequence 
stars and giants), eclipsing objects of both planetary and 
stellar nature (including white dwarfs), and configura- 
tions consisting of chance alignments with a foreground 
or background contaminant, as well as hierarchical triple 
systems. While these represent the most common and 
obvious configurations one can imagine, in principle there 
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Fig. 16. — Blend model of a background eclipsing binary with 
twice the nominal period of KOI-377.03. Top: Kepler observa- 
tions and best blend fit corresponding to two nearly equal mid G 
dwarfs eclipsing each other, and located 20.4 kpc behind the pri- 
mary. Bottom: Binned observations compared against the blend 
model in the top panel. The fit corresponding to a planet model is 
shown for reference. 



could also be more contrived scenarios that we have not 
thought of. These should be intrinsically much less likely, 
a priori, but can nevertheless not be completely ruled 
out. Therefore, we proceed below on the assumption 
that any such situations we have not considered have a 
small rate of occurrence, at least compared to the ones 
we have discussed explicitly. 

In order to provide the basis for an estimate of the 
confidence level for the planetary status of KOI-377.03, 
we describe here the calculation of the likelihood that 
the signal is due to a background blend involving either 
a stellar tertiary or a planetary tertiary, taking into ac- 
count the constraints on brightness and other properties 
indicated in the previous section. Because this type of 
calculation is likely to be relevant for other Kepler can- 
didates, we describe it here in some detail. 

The frequency of stars in the mass range permitted 
by BLENDER was esti mated using the B esancon Galactic 
structure models of iRobin et al.l (|2003l) . specifically for 
the R band, which is the closest available to the Ke- 
pler passband. We used an aperture of 1 square de- 
gree centered on Kepler-9, and we performed the stel- 
lar density calculations in half-magnitude bins of appar- 
ent brightness, accounting for interstellar extinction as 
we did in the BLENDER simulations, with a coefficient of 
0.5 mag kpc -1 in V. The range of allowed secondary 
masses for each magnitude bin was taken directly from 
Figure H2l for blends with stellar tertiaries, and from Fig- 
ure [13] for blends with planetary tertiaries. 

Using the density of stars in each magnitude bin, we 
calculated the fraction that would remain undetected 
after our high-resolution imaging (speckle and AO ob- 
servations), spectroscopy, and centroid motion analyses. 
The results are listed in Table [H The first two columns 
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TABLE 6 

Summary of blend configurations tested for KOI-377.03. 



False positive configuration a Result Blends ruled out 

Hierarchical triple with stellar tertiary, MS poor fits Yes 

• White dwarf tertiaries poor fits Yes 

Hierarchical triple with planetary tertiary, MS twin star Yes (imaging/spec. /centr.) 

Background EB with stellar tertiary, MS 

• With and without extinction good fits Not all 

• Giant secondaries P too short Yes 

• White dwarf tertiaries poor fits Yes 

• Twice the period poor fits/sec. eel. Yes 

Background EB with planetary tertiary, MS 

• Jupiters, Neptunes, super-Earths, with extinction .... good fits super-Earths 

• Giant secondaries P too short Yes 

• White dwarfs, with extinction good fits Yes (imaging/spec. /centr.) 

Foreground EB with stellar tertiary, MS poor fits Yes 

Foreground EB with planetary tertiary, MS poor fits Yes 

a 3 Gyr isochrone and solar mctallicity assumed for background and foreground stars. Orbits arc circular. Abbrevi- 
ations: MS — main sequence secondary; imaging/spec. /centr. — high-resolution imaging, spectroscopy, and ccntroid 
analysis; sec. eel. — secondary eclipses predicted but not observed; EB — eclipsing binary. 

give the Kp magnitude range of each bin and the mag- amples include contact binaries and ellipsoidal variables, 

nitude difference AKp compared to the target, calcu- in which the brightness changes continuously throughout 

lated at the upper edge of the magnitude bins. For the cycle, rather than presenting sharp transit-like events 

convenience the calculations for blends with stellar ter- such as we observe. Additionally, semi-detached systems 

tiaries and planetary tertiaries are listed separately. For would have eclipses that are too long and also of the 

the stellar tertiary case, column 3 reports the density wrong shape. We therefore exclude these from the tally, 

of stars obtained from the Besancon models, restricted With this adjustment, the frequency of eclipsing binaries 

to the mass range allowed by BLENDER as shown in Fig- capable of producing blends is 0.53%. Multiplying the 

ure [T2l Column 4 lists the maximum angular separation star counts in column 5 by this frequency, we obtain the 

Pmax at which stars in the corresponding magnitude bin total number of blends expected in each magnitude bin, 

would go undetected in our imaging observations, read which is reported in column 6. 

off from Figure [21 and taken at the center of each mag- Columns 7-9 are similar to columns 3-5, but for blends 
nitude bin. Centroid motion analysis rules out eclipsing involving star-planet pairs. Note that the range of al- 
binaries beyond 0'.'74, so p max is constant at that value lowed magnitudes is different in this case, as is the 
for the last few bins in which this provides a stronger range of secondary masses used to compute the densi- 
constraint than the imaging limits. The total number of tics in column 7 ( see Figure 1 131) . Follo wing the size 
stars of the appropriate mass range in a circle of radius ranges adopted by iBorucki et al.1 (|2010b| ). we consider 
Pmax around Kepler-9 is then given for each bin in col- three categories of transiting planets as potential com- 
umn 5, in units of 10 6 . We note that the secondary mass panions: super-Earths (1.25-2 i?©), Neptune-size planets 
and angular separation constraints together reduce the (2-6 i?©), and Jupiter-size planets (6-22 i?©, or equiva- 
number of background stars to be considered as poten- lently ~0. 5-2.0 i?j up ). Planetary tertiaries smaller than 
tial contaminants by a factor of 97,500 compared to the 0.3i?j up = 3.4 i?© are ruled out by BLENDER, as false pos- 
number that would otherwise be expected to fall within itives with such tertiaries can only reproduce the light 
the photometric aperture of Kepler. This is already in- curve if the secondaries are relatively bright, and those 
dicative of a significantly reduced chance of having a false would have been detected spectroscopically. This ef- 
positive. fectively eliminates all super-Earths, and a fraction of 
The intrinsic frequency of eclipsing binaries in the field the Neptunes. For the intrinsic frequency of transiting 
is a key ingredient in t he calculation, an d for this we have planets we have relied on the results fro m the first 43 
relied on the results of lPrsa et al.l (|2010t) , which are based days of Kepler observations as reported bv lBorucki et al.l 
on the Kepler observations themselves. These authors (|201 0bl. That census is unlikely to have missed many 
found the average occurrence rate of eclipsing binaries Neptune- or Jupiter-size planets (except ones with very 
among the Kepler targets down to Kp « 16 to be ap- long periods), although it may include false positives, 
proximately 1.2% across the entire field. This may be a so we consider the count to be conservative. Those au- 
slight overestimate because it counts as eclipsing binaries thors presented a list of 306 targets with at least one 
targets that are actually blended with a background bi- transiting planet candidate, and described another 400 
nary that is not a target, though the effect is likely small. targets with one or more transit-like signals that have 
There is little information available on fainter eclips- not yet been released. Of the 306 targets, 24% would 
ing binaries, so we assume here that a similar frequency correspond to Jupiter-size planets, and 23% to Neptune- 
holds. More importantly, many of these eclipsing bina- size planets. One may reasonably assume that the 400 
ries cannot produce signals such as that of KOI-377.03 sequestered targets contain a larger fraction of smaller 
because their light curves have the wrong shape. Ex- Earth- or super-Earth-size planets (which cannot mimic 
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the light curve of KOI-377.03; see Sect. 13.51) . for at least 
two reasons. Smaller planets are the main focus of the 
Kepler Mission, and they require more intensive follow- 
up efforts for validation, which is why those targets have 
not yet been made public. Secondly, these 400 targets 
are brighter, which makes smaller planets around them 
easier to detect. Consequently, the assumption of simi- 
lar Jupiter- and Neptune-size planet frequencies as given 
above, but for the entire sample of 306+400 targets show- 
ing transit-like features, is a conservative one. Scaling to 
the ful l sample of 156,097 Kepler targets (jBorucki et al.l 
l2010bf ) , we find an upper limit to the frequency of transit- 
ing Jupiter-size and Neptune-size planets of 0.11% and 
0.10%, respectively. 27 With these adopted frequencies, 
the resulting numbers of blends involving these types of 
planets are listed in columns 10 and 11 of Table [5] 

The total blend frequencies in each of columns 6, 10, 
and 11 are calculated as the sum of the individual fre- 
quencies in each magnitude bin, and the three columns 
are then combined in the bottom section of the table to 
yield an overall blend frequency (BF) of ~ 1.0 x 10~ 7 . 

This very small figure corresponds to the number of 
false positives we expect to find a priori for Kepler-9. 
However, we point out that this does not translate di- 
rectly into a false alarm probability, or equivalently into 
a confidence level that the candidate is orbited by a true 
super-Earth-size planet, as that requires knowledge of 
the rate of occurrence of such planets. For a random can- 
didate star in the field the rate of false positives relative 
to the rate of true planets (false alarm rate, FAR) can 
be written quite generally as FAR = N-pp/(Npp + N p ), 
where iVpp is the number of false positives and N p is the 
number of planets in the sample. Thus, the larger the 
number of planets we expect, the smaller the FAR. 

We consider Kepler-9 to be fairly representative of a 
typical target in the field in terms of its spectral type (so- 
lar), brightness, and background stellar density (a func- 
tion of Galactic latitude). In that case, the total number 
of blends can be taken to be approximately the product of 
BF and the size of the sample, or App = BF x 156,097 = 
0.016. The number of small planets expected in the sam- 
ple is of course not known, and determining it is precisely 
one of the goals of the Kepler Mission. 

If we accept a confidence level of 99.73% (3a) as being 
sufficient for validation of a transiting planet candidate 
(corresponding to FAR = 2.7 x 10~ 3 ), then the minimum 
number of super-Earth-size planets N p required in order 
to be able to claim this level of confidence is six. Even 
though N p is unknown, it is possible to make educated 
guesses as to what the minimum value would be in sev- 
eral ways, drawing on both theoretical and observational 
considerations. 

Ground-based Doppler surveys continue to push to- 
ward the detection of smaller and smaller planetary sig- 

27 As a check we may compare the above frequency of transiting 
Jupiter-size planets against results from the statistical study by 
IFressin et al. (200£j). These authors combined the findings of var- 
ious radial-velocity searches and folded-in the detections of very 
short-period systems detected by ground-based transit surveys, 
which are typically less common in Doppler searches. They found 
that 0.074% of solar-type stars have a transiting Jupiter-size planet. 
Our upper limit of 0.11% is consistent with this, given the presum- 
ably higher completeness of Kepler and the fact that a fraction of 
the Kepler candidates that is yet to be determined may turn out 
to be false positives once the follow-up is completed. 



nals. lLovis et al.l (|2008f ) have reported preliminary re- 
sults from a sample of some 400 FGK stars observed 
with the HARPS instrumen t on the ESO 3.6 m telescope 
(see also lMavor et a l. 2009), suggesting that as many as 
^30% of the targets may be orbited by close-in super- 
Earth- and Neptune- mass companions (5-30 M®) with 
periods up to 50 days. The peak in the period distri- 
bution seems to be around 10 days. By mak ing use of 
the CoRoTlux tool (jFressin et alj|2007l [2009h we simu- 
lated a sample of 156,097 stars in the Kepler field based 
on the Besa ncpn models emplo yed earlier, and used the 
results from lLovis et all (|2008|) to assign planets at ran- 
dom to each star, with a log-normal distribution of pe- 
riods centered at 10 days. Approximate planetary radii 
w ere inferred from the m asses using the structure models 
of iValencia et al] ()2007l ) . by drawing masses at random 
and assigning radii over the full range of compositions 
allowed by these models. We then retained only those 
in the super-Earth category, with R p < 2 R®. The num- 
ber of these planets that undergo transits in the Kepler 
field is calculated to be N p ~ 200. If we were to accept 
this estimate, the corresponding false alarm rate would 
be FAR = 8 x 10~ 5 . 

A similar Doppler survey of 166 G and K stars with 
the HIRES in s trum ent on the 10-m Keck I telescope 
((Howard et al] I2010D has provided additional insights 
into the rate of occurrence of small-mass planets. The re- 
sults suggest that approximately 18% of solar-type stars 
harbor planets in the range of 3 to 30 Mg with peri- 
ods under 50 days. A calculation analogous to that car- 
ried out above for the transit probabilities and conversion 
from planetary masses to planetary radii leads to an es- 
timate of N p w 120 for the Kepler field. If we adopted 
this lower estimate, the corresponding false alarm rate 
would be FAR = 1.3 x 10~ 4 . 

Pop ulati on synthesis studi e s such as those of llda &: Lin! 
(|2004T) and iMordasini et al.l (|2009D based on the forma- 
tion of planets by the core accretion process and subse- 
quent migration have provided tentative predictions of 
the properties of planets, including distributions of their 
masses, periods, and other characteristics. These the- 
oretical models seem to point to a sizable population 
of super-Earths that may be several times larger than 
the number of Neptune- or Jupiter-mass planets in those 
simulations. Scaled to the size of the Kepler sample, this 
would imply there could be several hundred small tran- 
siting planets. However, the authors caution that those 
results should be considered with great care as some of 
the physical ingredients in these models are still very un- 
certain. 

A further estimate may be obtained from the pre- 
liminar y Kepler results as reported by iBorucki et al] 
(2010b). Among the 306 targets listed there showing 
one or more periodic transit-like signals, 27 fall in the 
category of super-Earths. While it is true that these 
candidates have not yet been followed up and validated, 
and may therefore include some fraction of false positives, 
additional super-Earth-size planets are to be expected in 
the list of 400 unreleased candidates, which could make 
Np considerably larger. This is particularly true since 
the targets in the latter list are all brighter than those in 
the publicly available set, and therefore the proportion 
of small planets is likely to be higher because the transit 
events are easier to detect. Nevertheless, if we were to 
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accept that N p is as small as 27, then the corresponding 
false alarm rate would be FAR = 6 x 1CP 4 . 

There are caveats associated with each of the observa- 
tional estimates mentioned above that should be kept in 
mind. The Kepler results invoked in the previous para- 
graph are still preliminary, and although we regard our 
use of them to be conservative for the reasons described 
earlier, the true fraction of false positives in the Kepler 
sample remains unknown until all candidates have been 
followed up. The Doppler results are also preliminary to 
some degree, and are based on somewhat limited sam- 
ples of stars. It is also possible that a fraction of those 
Doppler candidates may turn out to be false positives, 
or given the sini ambiguity inherent in that technique, 
that some of them may have actual masses above the 
range considered for super-Earths. Additionally, there 
are uncertainties associated with the conversion we have 
applied between planetary masses and planetary radii, 
using theoretical models. Those uncertainties are diffi- 
cult to quantify given our present state of knowledge. 

For these reasons, added to the fact that despite our 
best efforts to assess the blend frequency there could still 
be some exotic blend scenario that we have overlooked, 
it is not possible to present a more definitive value of the 
false alarm rate. Nevertheless, the above estimates of the 
FAR based on consideration of all the blend scenarios 
that seem plausible to us are all sufficiently small that 
they give us very high confidence that KOI-377.03 is not 
a false positive, and they therefore validate it as a signal 
of planetary origin. We designate this planet Kepler-9 d. 

4. DISCUSSION 

Calculating the false alarm rate for targets with small 
signals such as Kepler-9 d is non-trivial because it de- 
pends crucially on the frequency of small transiting plan- 
ets, which may only be fully known at the conclusion of 
the Kepler Mission. Of the arguments for the expected 
value of N p presented in the previous section, the one 
that relies on the preliminary Kepler results themselves 
is the most conservative, and already makes it highly un- 
likely that we are in the presence of a false positive such 
as those explored in this paper. Furthermore, that esti- 
mate is based on results from only the first ^43 days of 
operation of the spacecraft. Continued observations over 
the next two years will surely increase the number of can- 
didates, which can only result in a larger confidence that 
the 1.6-day photometric signal is due to a true planet. 
Thus, we find the overall evidence for the planet inter- 
pretation very compelling. It is also worth noting that 
the Doppler surveys have found that a very large fraction 
of the smallest-mass planetary compa nions (as many as 
80%) are in multi- planet systems (e.g.. lLovis et alJ[2008t 
iMavor et al.l [200 9). Because of the presence of Kepler- 
9 b and c, this makes it considerably more likely that 
Kepler-9 d is also a planet than if it were the only signal 
in the system. Furthermore, one may expect a priori that 
a planet interior to Kepler-9 b and c would have a high 
probability of presenting transits. Indeed, with the rea- 
sonable assumption that the orbit of the inner planet is 
more or less coplanar with the outer two, the geometric 
probability of a transit at a period of 1.6 days would be 
close to 100%, instead of ~18% for random inclinations. 

The light-curve parameters we obtain for Kepler-9 d 
by modeling the photometry using the formalism of 



TABLE 7 
Derived properties of Kepler-9 d. 



Parameter Value a 

Orbital period (days) 1.592851 ± 0.000045 

Mid-transit epoch (BJD) 2,455,015.0943±°'oog3 

Orbital semimajor axis (AU) 0.02730±o;o°o43 

Transit duration (hours) b 1.971q^ 

R P /R* o.oi47tg;3SS 

Rp(r@) i-MlS:i" 

« 5-54±g;|J 

Equilibrium temperature (K) d 2026 ± 60 

a Values and uncertainties correspond to the mode and 1-cr confi- 
dence levels derived from the mode of the a posteriori distributions 
generated with the Markov Chain Monte Carlo algorithm. 

Defined here as the time interval between the first and last contacts. 
c The calculation of the normalized semimajor axis assumes the orbit 
is circular. 

^ Zero-albedo equilibrium temperature ignoring the energy redistri- 
bution factor. 

iMandel fc Ago] (|2002[ ) are summarized i n Table [71 and 
supers ede the preliminary estimates of iHolman et al.l 
(2010). The values and uncertainties were determined 
using a Markov Chain Monte Carlo technique, with four 
chains of length 10 6 each. Our fits used t h e non -linear 
fourth-order limb-darkening law of IClaretl ([2000). with 
coefficients for the Kepler band taken from the calcu- 
lations by A. Prsa referenced in footnote [231 For an 
adopted stellar radius for the primary star of R* = 
1.02 ± 0.05 Rq, the estimated size of this planet is R p — 
1-64!q;i4 i?® , which is among the smallest yet reported. 28 
The uncertainly is currently dominated by the photomet- 
ric errors, rather than the stellar parameters, and should 
improve as more measurements are gathered. The im- 
pact parameter of Kepler-9 d is very poorly constrained 
by current data, but should also become better deter- 
mined in the future. A mass determination for this ob- 
ject has not been made, and will be challenging given 
the small amplitude expected for the reflex motion of 
the star. The radial-velocity semi-amplitude would be 
only about 2.3 m s _1 assuming a similar mean density 
as the Earth. Velocity measurements are further com- 
plicated by the presence of the other two planets in this 
system. 

Many of the most interesting candidates found by 
Kepler will correspond to Earth-size planets, some of 
which are expected be in the habitable zone of their host 
star. For solar-type stars this implies reflex motions with 
radial-velocity amplitudes below current detection lim- 
its, making the spectroscopic measurement of the mass 
impossible. Other reasons may also hinder this type of 
validation, such as rapid rotation, chromospheric activ- 
ity, or even the faintness of the star. 

The information contained in the Kepler light curves 
on the shape of a transit-like event is a valuable asset 
for constraining the vast range of possible astrophysical 
false positives that might be masquerading as a planet. 
Here we have shown how modeling the light curve of a 

28 The size of Kepler-9 d is not significantly different from that of 
CoRo T-7b, which is R p = 1.68 ± 0.09 -R® according to Lcecr et al. 

(2009) , and was revised to Rp = 1.58 ± 0.10 by IBruntt et al.l 

(2010) . The measured radius of t he next smallest known p lanet, 
GJ 1214 b, is R p = 2.68 ± 0.13 (fCharbonnc au et alj|200Sh . 
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candidate directly as a false positive allows to rule out 
a significant fraction of blends involving background or 
foreground eclipsing binaries or star-planet pairs, as well 
as hierarchical triple systems, each of these possibly at- 
tenuated by the light of additional stars. The combina- 
tion of BLENDER with follow-up observations consisting 
of high-resolution imaging, spectroscopy, and an analy- 
sis of the centroid motion of the target leaves no room 
for a false positive in the case of Kepler-9b and c. These 
signals were previously known to correspond to bona- 
fide Sa turn-size planets bec ause they display correlated 
TTVs (jHolman et al.ll2010f ). Nevertheless, the exercise 
serves to show that, lacking that evidence, it would still 
be possible to validate them using the same techniques, 
thus supporting the general approach and justifying the 
application to the more interesting Kepler-9 d signal. 

Among the advantages of BLENDER for Kepler is the 
ability to predict the brightness and overall color of a 
blend in many different passbands. Brightness informa- 
tion for virtually every Kepler target is available from 
the KIC in the Sloan griz and 2MASS JHK S bands, 
as well as in the custom D51 passband (518.7 nm). Wc 
have shown earlier how these can be used to rule out cer- 
tain blend scenarios that might otherwise be viable. The 
detailed fitting of the photometry with a false-positive 
model provides additional discriminating power. Rough 
estimates of the properties of a background eclipsing bi- 
nary that can mimic a blend have sometimes been made 
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in previous transit surveys based simply on the apparent 
brightness of the object and a representative depth for 
its undiluted eclipses (such as 50%). While such config- 
urations may well reproduce the observed amplitude of 
a candidate light curve, not much can be said about the 
expected shape, which may be completely wrong. An ex- 
ample of this is seen for Kepler-9 d in Figure [TBI in which 
a back-of-the-envelope calculation of the depth one might 
predict from the brightness and ~50% deep eclipses of 
this binary may not be far off, but the detailed shape 
is not a good match to the observations, and BLENDER 
easily rules out this scenario giving a poor x 2 for the fit. 
Without these additional constraints on blend properties 
provided by the detailed light-curve fitting, the space of 
parameters open to false positives would be significantly 
larger, and more difficult to exclude by other means. 
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TABLE 8 

Blend frequency estimate for KOI-377.03 based on stellar densities and frequencies of eclipsing binaries and transiting planets. 



Blends involving stellar tertiarics Blends involving planetary tertiaries 



Kp range AKp Stellar density p max Stars EBs Stellar density p max Stars Transiting Jupiters Transiting Neptunes 

(mag) (mag) per sq. dcg (") (xlO 6 ) / EB = 0.53% per sq. deg (") (xlO 6 ) 6-15 R®, /j™ = 0.11% 3.4-6 R®, f Ncp = 0.10% 

(xlO 6 ) (xlO 6 ) (xlO* 5 ) 

(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) (11) 
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0.0359 


2674 
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Blend frequency (BF) = (0.0359 + 0.0339 + 0.0308) X 10~ 6 = 1.006 X 10~ 7 



